Casimir interaction from magnetically coupled eddy currents 
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We study the quantum and thermal fluctuations of eddy (Foucault) currents in thick metallic 
plates. A Casimir interaction between two plates arises from the coupling via quasi-static magnetic 
fields. As a function of distance, the relevant eddy current modes cross over from a quantum 
to a thermal regime. These modes alone reproduce previously discussed thermal anomalies of the 
electromagnetic Casimir interaction between good conductors. In particular, they provide a physical 
picture for the Casimir entropy whose nonzero value at zero temperature arises from a correlated, 
glassy state. 
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Spatially diffusive transport is a basic physical phe- 
nomenon that has been studied with a wealth of meth- 
ods. For example, the equation for heat conduction was 
solved by J. Fourier by his transformation that provided 
later the framework for quantum field theory. Remark- 
ably, the Fourier modes of the diffusion equation itself, 



d t A = DV 2 A, 



(1) 



do not fit into a simple quantum field theory because they 
have purely imaginary frequencies. The quantization of 
(over)damped modes can be done, however, with an al- 
ternative approach, the 'system+bath' paradigm of dis- 
sipative quantum dynamics. In this picture, the obscrv- 
ables relevant to a mode (specified, e.g., in Fourier space 
by its wavevector q) are damped in time because they 
strongly couple to a system with infinitely many degrees 
of freedom that are not directly accessible ('bath') [J]. 
The additional fluctuations the bath couples into the sys- 
tem, compensate at equilibrium for dissipative loss [2], 
even at zero temperature. This implies that a quantity 
like the zero-point energy of a damped mode is no longer 
given by the usual ^hujq and must be redefined and rein- 
terpreted. Indeed, the ground state of the combined sys- 
tem+bath is, in general, entangled: the corresponding in- 
teraction Hamiltonian is responsible for the change in the 
zero-point energy relative to the decoupled systemj^, 
Clearly this has an impact on all phenomena connected 
with fluctuation energy, for which the paradigmatic ex- 
ample is the Casimir effect 0. 

In this letter we discuss a specific example of 'diffu- 
sive modes': eddy (Foucault) currents in two metallic 
plates Q. We address their quantum and thermal fluc- 
tuations and their role in the electromagnetic Casimir 
interaction Similar to the approach of Ref.0], we iso- 
late the contribution of eddy currents among all other 
modes. We show that these modes quantitatively repro- 
duce the "unusual features" of the Casimir force between 
good conductors at finite temperature (as predicted by 
Lifshitz theory), that have been intensely debated for 
several years [a, |9( . 

To begin with, consider a metallic bulk described by 



a local dielectric function in Drude form: e(u>) = 1 — 
fl 2 /[uj(uj + ij)] (0: plasma frequency, 7: damping rate). 
This system allows for a continuum of damped, charge- 
less modes with a (transverse) current density and (dom- 
inantly) magnetic fields at frequencies of order 7 or be- 
low: they are called 'eddy' or 'Foucault currents' and are 
at the base of phenomena like magnetic braking or the 
familiar induction ovenQ. The 'bath' for these modes 
is provided by the phonons in the metal or by impurity 
scattering. On time scales slower than I/7, the vector 
potential satisfies Eq.([TJ) with an electromagnetic diffu- 
sion constant D — 7A 2 (A = c/fi: plasma penetration 
depth, A w 20 nm for gold) . A mode analysis in Fourier 
space yields imaginary frequencies 07 q = — i£ q « —\Dq 2 
for q <C I /A: due to the coupling to the electromag- 
netic field, the damping of eddy currents becomes weaker 
(compared to 7) and dispersive. 

Let us now split the metal in two half-spaces, inserting 
a vacuum layer of thickness L ( "cavity" ) . While the eddy 
currents cannot leave the medium, the associated elec- 
tromagnetic fields do; they cross the cavity in the form 
of evanescent waves. The corresponding magnetic fields 
provide Johnson noise that dominates, at low frequencies, 
over the fluctuations of other modes confined within the 
cavity [Tol II, [T3|. The characteristic frequency is the 



so-called Thouless energy jl 31] 

h£ L ee HD/L 2 



(2) 



and the scattering rate 7 provides a strict upper limit to 
the £ q . Typical numbers for gold at room temperature 
are h( L ~ 20 K and hj ~ 500 K at L = 100 nm. We 
focus in the following on the polarization with the elec- 
tric field parallel to the metal surface (TE). Indeed, the 
orthogonal case (TM) is associated with a surface charge 
that effectively decouples medium and vacuum. 

For an overdamped mode, the 'system+bath' approach 
predicts a zero-point energy^, [Til 



(3) 



where A is a cutoff frequency related to the bath response 
time (HE!. The corresponding Casimir energy captures 
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the shift in the (imaginary) mode frequencies £ g = £, q (L) 
due to the magnetic cavity field correlating the eddy cur- 
rents in the two plates. We therefore sum this expression 
over the eddy current mode continuum and subtract the 
reference system of two isolated half-spaces The dif- 
ference of mode density per area, p(£; L) is related to 
scattering phases at the metal- vacuum interface (see Ap- 
pendix) : 



d 2 k 

7r(27r) : 



Im log [1 - r£(-i£ + 0)e~ 2Ki ] 



where k = y/k 2 + (£/c) 2 and rk the TE reflection co- 
efficient. This shows power laws in £ below and above 
the Thouless frequency £l, providing the asymptotic be- 
haviours [FigQ] 



E D (L) 



const. - hjLX 3 log(A/7), L < A, 
n{t L c/L) l ' 2 L- 2 log(AL/c), L » c/r 



(5) 



The associated pressure —dEjj/dL shows a repulsive 
Casimir effect flOl. 11 . ll2T |. provided the cutoff is suffi- 
ciently large (a few times 7). The repulsive character 
can be understood from the Lenz rule for current-current 
(diamagnetic) interactions @ . Eq.© does not show the 
scaling hc/L 3 characteristic of perfect reflectors [H and 
is much weaker because of the finite bandwidth < 7 of 
the relevant mode continuum. The pressure vanishes for 
7^0 because dynamical fluctuations induced by the 
bath are suppressed. At small distances, En(L) essen- 
tially scales with the 'missing volume' L x (plate area) 
and the bulk mode density ~ 'y 1 / 2 £)~ 3 / 2 (at £ ~ 7) of 
the eddy current continuum. Their zero-point fluctua- 
tions exert, in this regime, a constant kinetic and mag- 
netic pressure. 
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FIG. 1: Casimir pressure (T = 0) between two metallic plates: 
contribution of eddy currents alone and of all electromagnetic 
modes (Drude metal) in the TE polarization. Normalization 
to the perfect reflector, 7r 2 7ic/(480L 4 ), distance L in units of 
the the plasma wavelength (A p = 2nX). Dashed lines: asymp- 
totics Eq.©. We take 7 = 0.08 SI and a cutoff A = 57. 

Let us now raise the temperature. Eddy currents are 
thermally excited at quite low temperatures because of 
their low characteristic frequencies: for good conductors 
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FIG. 2: Eddy current mode density p(u>; L) for two plates at 



(4) distance L (7 = 10~' 5 SI, 



A p = 2tv\, 



a'/L 2 = £l\p(0;L)\ 



Eq.©]). Frequency normalized to the Thouless frequency 
£l = D/L 2 . Thick gray line: mode density used in the Lif- 
shitz theory with the Drude dielectric function. TE polariza- 
tion only. 



at room temperature we have H£l <C k B T < 717 already 
at L = 100 nm. We find indeed that they give a signif- 
icant contribution to the Casimir interaction. Consider 
first the Casimir entropy. The contribution due to eddy 
currents can be calculated in terms of a real frequency 
mode density ('DOS') p(u;L) [FigJ5] using the following 
Kramers-Kronig like relation: 



p(uj;L) 



TT ^ 



(6) 



Fig|5] shows that the DOS is significant in the frequency 
range to ~ . . . 7. One also sees that in this range, 
the DOS of the electromagnetic Casimir effect between 
Drude metals (gray line) is entirely given by Foucault 
currents. The zero frequency limit is 



p(LU^Q;L) 



2 In 2 - 1 
8tt 2 D 



< 0. 



(7) 



Given this DOS, we get the free energy J-(T, L) by mul- 
tiplying with k B T\og[2smh(huj/(2k B T))] and integrat- 
ing over oj (see Appendix). The entropy, S(T,L) = 
—dT/dT, is cutoff independent and becomes [see also 
FigE] 



TT 2 k B T 



p(0;L), 



S_ = I 3 h 

k B 1 C(3)/(£/A) 
16ttL 2 



k B T^Tii Ll 



(8) 



Soo(L), k B T^>n^ L , 



where the dimensionless function f(L/X) — > 1 for L ^> 
A. The characteristic temperature scale h^/k B has 
been noticed previously from the Lifshitz theory with a 
temperature-independent scattering rate [1, [l(| ■ A neg- 
ative entropy is not surprising here since we are consid- 
ering a difference, with the entropy for infinitely sepa- 
rated plates subtracted. The sign means that one plate 
is acquiring information (entanglement) about the con- 
figuration in the other one through the electromagnetic 
interaction. The entropy vanishes linearly as T — > 0, in 
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agreement with the Nernst heat theorem (third law of 
thermodynamics). In this limit, modes with frequencies 
above k B T/h 'freeze' to their ground state, do not con- 
tribute to the entropy, and a unique ground state for the 
system remains at T = 0. We have checked that the 
free energies at low T due to eddy currents alone and in 
the full Lifshitz theory [3, G3] coincide in their first two 
terms (T 2 and T 5 / 2 ). 
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FIG. 3: Eddy current Casimir entropy vs temperature T. 
Dashed: Eq.([SJ) for T — > 0. We normalize to the perfect reflec- 
tor entropy at high temperatures, a/L 2 = (^(3)fcs/(167rL 2 ). 
Arrows: crossover temperature ksT — h^L = TiD/L 2 . 

Another scenario emerges with a T-dependent scatter- 
ing rate when the ratio £l(T)/T vanishes as T — > 0. 
This happens for example in a 'perfect crystal' fl8j where 
j(T) ~ T n (n > 2). In this case the Foucault modes 
stay in the high-temperature regime and give a univer- 
sal contribution S^L) ~ —f(L/\)/L 2 to the entropy. 
The same 'entropy defect' has been found in previous 
work on the Lifshitz theory with the Drude modelfLsj 
and traced back to evanescent TE modes (T^]. Neverthe- 
less this does not conflict with thermodynamics because, 
in this limit, the ground state becomes infinitely degener- 
ate: for D -> 0, the DOS p(0; L) -» oo in Eq.flf}, and the 
spectrum collapses to zero frequency because the char- 
acteristic frequency scale £z, — > 0. The diffusive modes 
then become static current loops interlocking with mag- 
netic field lines that permeate the metallic bulk, simi- 
lar to the frozen magnetization in an ideally conducting 
medium. This is a glassy state with a nonzero entropy 
('Foucault glass'). As a check of this idea we have for- 
mulated a Hamiltonian field theory for a glass of static 
currents and have calculated the change in entropy per 
area for two magnetically coupled half-spaces, recovering 
a result that perfectly coincides with Soc(L). It is clear, 
of course, that the Foucault glass is the result of an over- 
idealized characterization ('perfect crystal'). Indeed, the 
Lifshitz entropy vanishes with T — > when additional, 
even small, relaxation mechanisms like Landau damping 
are included 0. 

Let us now consider the eddy current contribution to 
the Casimir pressure at T > 0. The debate around 
the electromagnetic Casimir interaction 0, 0] has re- 
vealed that the large-distance pressure between Drude 



metals is effectively halved compared to perfectly re- 
flecting mirrors, even for an infinitcsimally small dis- 
sipation rate. This is illustrated in Fig0] where the 
temperature-dependent part of the Casimir pressure is 
plotted for the Drude model (gray lines vanishing as 
L > At = hc/(2k B T)). The Drude model thus can- 
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FIG. 4: Thermal Casimir pressures (TE polarization 
only) vs. distance L, normalized to the perfect reflector 
[C(3)fc s r/(87rL 3 )]. In the Drude model (gray, 7 = 0.08 0), 
the pressure is repulsive and vanishes for L 3> Hc/(2kBT) 
(black arrows). Blue dots: Drude model with eddy current 
pressure subtracted. Black: plasma model. Gray arrows: 
electromagnetic diffusion length [hD / (kBT)] 1 ^ 2 . 

not describe the limit of a perfect reflector (dashed line 
in Fig (4|), even if we take 7 — > or A — > ('ideal con- 
ductor'). Models that go over to the perfect reflector (as 
L > At) are superconductors (20I. w\ and the so-called 
'plasma model' (lossless dielectric function with 7 = 
right from the start, solid line in Figf4]). The thermal 
pressure from eddy currents accounts nearly completely 
for the difference between the Drude and the previous 
models If], U , 13] : they thus represent very precisely the 
object under debate. In particular eddy currents provide 
a simple physical explanation why in the Drude model, 
the thermal Casimir pressure sets in at unusually short 
distances (compared to the thermal wavelength, arrows 
in Fig[4]): this happens indeed when the dissipative en- 
ergy scale crosses over from the quantum into the clas- 
sical regime, HD/L 2 ~ k B T The absence of eddy 
currents in the plasma model is actually not surprising 
since it behaves like a superconductor and expels low- 
frequency magnetic fields [22] • The previously discussed 
'perfect crystal' limit (D — ► of the Drude model), on the 
contrary, corresponds to the ideal conductor that does 
not show the Meifiner effect [22 1. 

Recent measurements of the Casimir interaction be- 
tween gold-coated bodies (a plate and a sphere) at room 
temperature have been argued to favor a theory within 
the (lossless) plasma model [23]], but these results do not 
appear to be commonly accepted yet. In view of our 
analysis, a similar prediction could equally emerge within 
the Drude model if the contribution of eddy currents 
were somehow reduced, keeping only the other (photonic) 
modes. How this should be implemented physically, re- 
mains to be understood in detail. Potential mechanisms 
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should take into account, however, that in other contexts, 
thermally excited eddy currents have been observed, in 
agreement with the usual thermodynamic theory 24 1 . 



To reduce the eddy current Casimir interaction, we 
have considered a calculation at constant entropy, noting 
that also the thermalization rate of low-frequency modes 
is very small. This approach leads, however, to a signif- 
icant overshooting of the Casimir pressure at large dis- 
tances (very low specific heat according to Eq.©). Any 
mechanism that gives a nonzero real part to the frequen- 
cies of eddy current modes (see e.g. Ref. would re- 
duce their contribution at low temperatures because they 
then 'freeze out' at fcgT <C Re foujk- This is not likely for 
current room temperature experiments, however, as esti- 
mates based on finite-size effects show. 

In summary, we have discussed an example of quantum 
thermodynamics for an overdamped field, calculating the 
Casimir energy due to (suitably generalized) zero-point 
energies and its dependence on temperature. This applies 
to a simple physical system of two metallic plates that has 
been controversially discussed for several years by now. 
We have considered diffusive (eddy or Foucault) currents 
in normally conducting plates and their correlations via 
magnetic fields that lead to a repulsive pressure. These 
modes have a range of characteristic frequencies that is 
very different from photon modes. They cross over from 
a quantum to a thermal regime in a range of parame- 
ters that is accessible with current experiments, namely 
when the energy scale of diffusive transport (Thouless en- 
ergy hD/L 2 ) exceeds fc^T. The finite-temperature quan- 
tum field theory of these modes gives significant entropic 
corrections to the Casimir pressure that set in at dis- 
tances shorter than the thermal photon wavelength and 
explain quantitatively the large thermal corrections to 
the Casimir pressure in the Lifshitz theory for Drudc met- 
als. We reproduce the nonzero Casimir entropy for this 
system and suggest a physical understanding of its un- 
usual behaviour. In the particular case that the damping 
rate becomes infinitesimally small at low T, an entropy 
defect remains that has been interpreted in terms of a 
glassy state of quasi-static Foucault currents, the Nernst 
theorem being not applicable. This limit does not recover 
the lossless case (plasma model) because there, magnetic 
fields are expelled from the bulk, as in the London de- 
scription of a superconductor. Similar conclusions could 
also apply to other classes of modes (like overdamped 
coupled surface plasmonsQ) providing insight into sim- 
ilar controversies [25| and perhaps a route to engineer 
the Casimir force, similar to the suggestion inQ- From a 
theoretical point of view, our analysis illustrates that the 
unusual behavior of quantum field theory for normally 
conducting mirrors finds a clear physical justification in 
the framework of open quantum systems, glassy dynam- 
ics, and superconductivity, and paves the way for further 
theoretical and experimental investigations. 
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Appendix: eddy current mode density 

For the density of modes at real frequencies, p(cu;L), 
of the electromagnetic Casimir effect we use the following 
prescription: the free energy (per unit area) is written in 
the form 



F 



du p(u;L)f(iv) 



(9) 



where f(uj) = kBT\og[2smh(huj /2ksT)] is the free en- 
ergy per mode. For two parallel, identical plates, this 
gives 



D„(k,w) = l-rjl(ft,w)« 



(11) 



The integration is over the wavevector k parallel to 
the plates, and the two polarizations p = TE, TM are 
summed over. The r^(k,oj) are the reflection coefficients 
for the electromagnetic field from a single half-space, and 
k is defined after Eq. (fT4")) below. Passivity does not allow 
for poles and zeros of the dispersion function D AI (k, uj) 
in the upper half of the complex plane. Its analytical 
continuation into the lower half plane does show, how- 
ever, poles, zeros and branch cut singularities: the zeros 
(poles) actually define the complex mode frequencies of 
the two mirror system (of two isolated mirrors), respec- 
tively. Branch cuts arise by taking a difference of mode 
densities in the continuous part of the spectrum. The 
expansion of the mode density p(u; L) over these singu- 
larities, displaying explicitly the contribution of a branch 
cut along the negative imaginary axis, reads 14 1 



p(uj;L) = 



7T LJ + i£ 



Pother modes (u;L) (12) 



where we have introduced the (real- valued) mode density 
along the branch cut p(£; L). The first term is Eq.(6) of 
the main paper. 

In order to calculate the mode density p(£; L), we spe- 
cialize to a specific form for the reflection coefficients: 



5 



namely the Fresnel form with the Drude model 

n 2 



1 



lo(lo + ry) 



(13) 



involving the plasma frequency f2 and 7 the scattering 
rate for the current density. In the TE and TM polariza- 
tion, respectively Q 



r T E{k,uj) 



e{m)n-n m 
r T M{k,uj ) = - r - — ■ 14) 



where k = y/k 2 — (lu/c) 2 and K m = \Jk 2 — [lu / c) 2 e{lu). 
The branch cut at negative imaginary frequencies lu = 
— i£ is associated with the square root K rn ; it connects 
the solutions of the equations 

k 2 - e{lu){lu/c) 2 = 0, -00 (15) 

The second equation easily gives w = — 17; the first equa- 
tion is a third order equation of which we pick the purely 
imaginary root lu = — i£k- (The other two correspond to 
the plasma edge above which electromagnetic waves pen- 
etrate into the bulk of the mirrors.) In the limit k <C f2/c, 
we get 



= -ifk ~ -iDk 2 



(16) 



where D = jc 2 /tt 2 is the diffusion coefficient. Both 
reflection coefficients and hence the dispersion func- 
tion (fTTj) show a branch cut for lu 6 [— i£k, _ 17]- The 
mode density /5(£; L) along this cut is found from the ar- 
gument principle, taking an integration contour around 
the negative imaginary axis. This gives 

p(i;L) = -d c Tm^2 /^2 Im IogD M (k,-i£ + 0) 

(17) 

where the +0 prescribes the side of the branch cut in the 
right half plane (setting the sign of the wave vector n m ). 

The same result can also be obtained using a scatter- 
ing approach. In this case, we begin with the dispersion 
relation in a bulk Drude metal, 







c 2 q 2 — lu 2 £(lu) 



c 2 q 2 



1 



O 2 



c 2 q 2 



(18) 



LU Z LU + 17 

where e(tu) is again the Drude model. Both dispersion 
functions in the numerator and denominator of this ex- 
pression are a simple re-writing of the wave equations 
after Fourier transformation. The poles of Eq. (fT8|) il- 
lustrate the subtraction operated in the bulk dispersion 
energy: the free space photon continuum and locally 
damped currents (not coupled to the long-range electro- 
magnetic field). Looking for the solutions to Eq. (fT5)) . 
three modes are found, of which one always occurs at 
imaginary frequencies cj q = — i£ q . Its limits for small 
and large q are 

-\Dq 2 



q <C VL/c : w q 
q > fl/c : w q 



(19) 
(20) 



where D is again the electromagnetic diffusion constant. 
In view of planar surfaces for the metal, we label the 
modes by the wavevector k projected onto the plane z = 
and the frequency — i£. The perpendicular wavevector 
k z is then fixed by the dispersion relation (jTSJ) - 

The electric current density associated with diffusive 
waves cannot escape into vacuum (except for the dis- 
placement current). Hence we get total internal reflec- 
tion with a coefficient of modulus unity. In the vacuum 
gap, the (transverse) electromagnetic field varies expo- 
nentially with a decay constant k 2 = k 2 + £ 2 /c 2 . The 
single-interface reflection coefficient is obtained applying 
the usual boundary conditions for the electromagnetic 
field. In the TE-polarization, we find 

r D (k,-iO = -^r = -rTE(k,-^ + 0) (21) 

K — lfc z 

With two interfaces at z = ±L/2, odd and even modes 
with respect to the mid-plane of the cavity vary as sinh kz 
and coshfcz for \z\ < L/2. In the medium, we find that 
the reflection coefficient rjj is modified by a phase shift 
c" 2 " 5 , and this gives access to the change in the mode 
density (Le., the DOS p(k, £)) of the diffusive wave con- 
tinuum [26j . The phase shifts are measured with respect 
to the configuration where the mirrors are placed an in- 
finite distance apart. Alternatively, we can compute the 
transmission coefficient across the vacuum gap and get 
the DOS from its phase[27j. Summing over even and 
odd modes we get a phase shift 



_ m = l-{r* D ) 2 e- 2 ^ 
l—r 2 D e~ 2KL 



(22) 



The mode density along imaginary frequencies, per k- 
vector and polarization, is thus 

p(k, = -i^Im log [1 - r 2 D (k, -iO e- 2Ki ] (23) 

Integrating over the parallel wavevector and summing 
over the polarizations, we recover Eq. (|17[) . 
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